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ABSTRACT 


Numerical weather analysis relies heavily on numerically computed 
winds both for direct output and internally for calculating advection and 
other purposes. It is the intent of this paper to compute, analyze, and 
compare the geostrophic, the non-divergent, the divergent, and a com- 
bination of the latter two. Where there occurs a significant difference 
between these numerically derived winds, an attempt is made to com- 
pare each to the observed winds. This study is made for all standard 
levels through 300 mb. Data and several computer programs were fur- 
nished by the U. S. Navy Fleet Numerical Weather Facility, Monterey, 
California. 

The writer wishes to express his appreciation to Professor George 
J. Haltiner of the U. S. Naval Postgraduate School for his guidance, 
contributions and encouragement in this work. 

Appreciation is expressed to Lieutenant Commander Mildred J. Seljos, 
United States Navy for aid and advice in programming. Further thanks 
are extended to Lt. R. F. Alden, U. S. Navy, and Lt. G. C. Rosenber- 
ger who provided programs for computation of horizontal velocity diver- 


gence and for advection of temperature. 
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Vg 
Ve 
VW 
Vx 
Vt 


Va 


I¢ 
DGD¢ 
IGI¢ 
DADV 
EADV 
TADV 


XADVT 


re 


u¢ 
v¢ 


D¢ 


finite difference Laplacian vA s(n, tir Ai, hie reve 
finite differnece a fai. at aa oe (pA (ps B?,) | 
geostrophic wind 

non-divergent wind 

Ve 

divergent wind 

total wind 

actual wind as measured from soundings 

dummy variable 

wind speed of V¢ 

difference in wind direction, Dg - D¢ 

difference in wind speed, Ig - I¢ 

advection of absolute geostrophic vorticity by Vg 
advection of absolute stream function vorticity by Ve 
advection of absolute stream function vorticity by Vt 
advection of temperature with Vx 

time 

u component of V¢ 

v component of V¢ 


wind direction of V¢ 





1. Introduction 

One of the primary products of any numerical weather facility is 
wind data for operational useand for further numerical computations 
such as advection. It is the purpose of this paper to analyze and 
compare several types of numerically derived winds. The types con- 
Sidered are the geostrophic, Vg, the non-divergent, Ve, the divergent, 
Vx, and finally the total wind, Vt, which is a combination of the di- 
vergent and non-divergent parts of the wind. These winds are cal- 
culated by solution of appropriate equations using data from the U.S. 
Naval Fleet Numerical Weather Facility, Monterey, California (FNWF). 
Computations were made on the U. S. Naval Postgraduate School CDC 
1604 digital computer. 

Comparison of these winds is made at all standard levels through 
the 300-mb level by calculation of speed and geographical direction on 
separate maps for each level covering the Northern Hemisphere. Dif- 
ference fields of direction and speed are made between each wind and 
the geostrophic wind. Vertical cross sections were computed and 
studied for each of the winds and for the difference fields. Further com- 
parisons and analyses are made by using these three winds to advect 
absolute vorticity. Advection of temperature was also computed using 


Vx. 





2. Background 


The geostrophic wind is computed from the well known formula 
My ESR E - Beg) 
jg f \ oy = (1) 


utilizing a standard program of FNWF. This program generates fields 
of geostrophic wind directions in degrees and speed in knots. All com- 
putations were done on the FNWF square grid centered on the North 
Pole and consisting of 3968 with a grid interval of 381 km at 60N. 
In this program the sine of Jatitude, as required in computation of the 
coriolis parameter, f, was set equal sine 15 degrees. The latter points 
were not utilized in the evaluation of results. 

Geostrophic wind, Vg, is utilized as a basisofcamparison for 
later derived non-divergent, Ve, and total, Vt, winds. 

The Helmholtz theorem allows expression of the wind as the sum 


of the divergent and non-divergent components. 


Winds of this form will be her *tesi referred to as total wind, Vt.. 
The potential function x for the divergent part of the wind is obtain- 


ed through relaxation of the Poisson-type continuity equation 


eats =. (3) 





The stream function, ¥, for the non-divergent portion of the wind is 


approximated by solution of the balance equation 


ty \2 = 2 
vr vrv eve -#l (Sea)- $8 i v9. (4) 


Here, @is the geopotential, 
The vertical component of vorticity, 0, is the Laplacian of the 


stream function 


Ceres (5) 


and aksolute vorticity, 7 , is the sum of fand the vertical component 


of the earths vorticity, fj 
i ee (6) 


Vertical velocity, w, is computed by the diagnostic equation 


PE ia Ae IE | 

Vout aed vi J (7) Joy T) a cere) (7) 
developed by Haltiner, Clarke, and Lawniczak Cue 

For convenience the stream function has been redefined so as to be 


in units of height, 


Y= ZTE (8) 


The geostrophic advection of absolute geostrophic vorticity is com- 


puted from the expression 


aes 
Te J (0, Ng) (9) 





and the advection of n, by 


ee) (10) 
= J(E, Ne). 


The advection of ~ by Vx is computed from the formula 


We V Ne = Ux alle + Vx = - 


(10) and (11) are mer to represent advection of nN. by Vt and 





compared to similar advection using Vg and Ve. 





3. Procedures 


The sequence of calculations programmed started with solution of 


the geostrophic wind equation in finite difference form 
Ginn, ee J j a | 
\V = aa) (D; - 0; Je - (Dis —Di-)4 . (12) 


Fields of Ug and Vg relative to map projection were generated and convert- 
ed into fields of wind direction in geographical headings at each grid 
point. Likewise, fields of isotachs were produced and printed in con= 
toured map form. These calculations and maps were made for all stand- 
ard levels to 300 mb for this and all subsequent programs. 

The balance equation, (4), as programmed by M. Reese of FNWF, 
was then solved and the fields of E values were stored on magnetic 
tape. The iterative scheme of this program does not always converge at 
300 mb or has an excessive running time, both of which limited the data 
produced. The E values were then used as input to the geostrophic 
wind program with the appropriate substitution of f for f. These "balanced" 
winds were computed and _ stored in the same manner as Vg. 


Divergence fields are taken as input for solution of (3). Here, as 


with Y, X is redefined as 


(3) 





so that the velocity potential has units of height and therefore is com- 
patable with E and D forus2i1 (12) with minor changes. Expressing 


the first term of equation (3) in finite difference form gives 


moe td = Qu 
Vx "am Op (14) 


which may be solved by use of an FNWF Poisson relaxation program. 
The derivitive ~3u/ 3p which is the horizontal divergence of the wind 
is evaluated by a finite difference approximation as programmed by 
Alden and Rosenberger [2] . 

Tht.” fields are then used to compute the divergent wind component 


by the finite difference equation 
x) as ~ Jti Y Je - | 
Ie ((Z},- Rik + (Ze (15) 


Again the u and v components are stored and maps of direction and of 
isotachs are produced. 
Now the components of Ve and Vx are added and a field of these com- 


bined components produced as follows 
Ugyt Ux = Ut (16.a) 


Vat Vx =Vt (16.b) 


As before these component fields are stored and maps of each level print- 
ed. 

Programs to aid comparison and analysis are now written. Analysis 
of directional differences is made by using the direction fields of Vg as 


a basis and subtracting directions of Ve and Vt from it. This proved to 








be a challenging problem for several reasons. First, directions are 
by nature discontinuous at singular points on isogons maps and also 
in general at 360 degrees. Second, direction tends to vary widely 

in the vicinity of small scale features that are not concentric between 
the various fields compared. 

The discontinuity problem was solved successfully with a rather 
intricate program. However, the singular points and non-concentric 
small scale features still caused a large gradient of direction and 
difficult contouring to plot or read. This problem was alleviated by re- 
placing a value of minus zero at all points on the directional difference 
maps at which the wind speed was less than five. knots. This improved 
appearance of the maps and facilitated evaluation. Moreover the min- 
us zero values still indicated those areas where singular points had 
occurred. The 5=-knot criterion proved to be of no detriment in subsequent 
use. 

Speed-difference fields, IGIE and IGIT, are produced by the same pro= 
gram as directional difference. The computation of speed differences in= 
cluded a statistical analysis and grid maps of output included the pillow, 
P, and root-mean-square difference, RMSD. These terms are defined in 
APPENDKX I. 

Another series of programs were written to calculate the advection of 
absolute vorticity, 7 , by the various winds previously computed and to 


compare the results. 





¢ 


x 


First, fields of D values are used as “input to evalluate ‘the geostro- 


phic vorticity yn 





(17) 


with an FNWF program. Then the reuters fields from (17) are used to 


Mt f=- Jat mr Teor ”, ) (18) 


which is plotted on standard grid maps. This program is then internal= 


obtain 


ly modified to solve 


Ne = te ee 


(19) 


and 
Ay — ioe J) (E He): (20) 


Next previously stored components of Vx are used to compute the 


advection by the divergent wind as follows: 


A, i} _ atm a Me Nin)t i m% (Mi th) | | (21) 


This result is added to (20) to produce fields of TADV. 

Next these three advection schemes are compared using DADV as a 
basis and forming difference fields similar to those used for direct com- 
parison of the winds. 

Typical computer operating time including efficient handling of tapes 
and tee of output on an IBM 1402 printer is shown in Table I. 


*2es shown are for one reporting period of one day at all levels. 





’ Table I 





Program Run time in minutes 

Solution of Balance Equation, (4) 46 

Wind speed and directions from various inputs 18 €ach) 

Solution of equation (14) for x 30 

Solution of (16.a) and (16.b) 12 

Comparison of winds 35 

Advection of absolute vorticity with all winds ops 

Comparison of advections 20 

Preparation of grid field for FNWF plotter 4 (each) 
Total 200 minutes. 


It is apparent that operating time is a significant factor in restricting the 


amount of data collected. 





4. Results 

Results and conclusions are based on analyses of data selected 
from the period February 10 to 20, 1964. Only the United States is 
used in figures as it represents an area Of reliable data. Usinga 
small area for figures also allows display of finer detail. For illustra- 
tion 0000Z 16 Feb 1964 was selected as it contains a wide range of 
synoptically interesting features over the United States. 

Figure 1 depicts the pressure analysis (dashed lines) superimposed 
on the stream function analysis at 500 mb. The relative height differ- 
ence at the semi-enclosed low center south of the Great Lakes is typ- 
ical of the 500-mb level throughout the period studied. Figure 2 shows 
stream function analysis of the same period. Comparison with pressure 
height analysis iidicates that the patterns are similarly related at other 
standard levels. 

Wind speeds computed from 500-mb height and stream function analy- 
sis are shown in figure 4. Dashed isotachs enclose centers of maximum 
and minimum speeds. Speeds for other standard levels are shown in 
figures 2 and 3 where dashed lines are isotachs. Large numerals indicate 
values of minimum and maximum speeds. Analysis of speeds of these 
two winds shows there is a significant difference between them. Dashed 
lines in figure 5 (lower) are isolines of difference in wind speeds, geo-: 
strophic minus non-divergent. Similar analysis of other levels are shown 
in figures 6 and 7. 


10 





Figure 5 (upper) shows differences in direction for 500 mb, Dg - De. 
Solid lines show the 500-mb pressure height analysis, alternate dashed- 
dot lines arecontours of negative 5-kt speeds, and dashed lines are 5-kt 
contours of positive difference. The area marked > 30 is almost com- 
pletely within an area of wind speeds of less than 5 kt. Difference in 
directions are not shown for other levels but are easily located by refer- 
ence to figures 2 and 3. Pillow and RMSD values for the hemispheric 
maps is given in figures 5 through 7. Note that P indicates that on the 
hemispheric scale Ve speed is lower at upper levels, e.g. by a factor 
of approximately 10 per cent at 300 mb. Except in regions of very sharp 
curvature or near singular points, DGDE, figure 5 (upper), is generally 
less thantendegrees at this and all other levels. A more detailed 
analysis was made of gridpoints and actual winds from 185 soundings 
in areas where DGDE was greater than 10 and both Ve and Vg exceed 
5 kt. The graph in figure 8 shows therealationship found. It must be 
pointed out that interpolation within grid points was required to use 
actual soundings and that actual! soundings are transmitted to the nearest 
10 degrees only. Inan effort to at least partially reduce effects of fric- 
tion and local terrain effects the surfice winds were not included in the 
data used for construction of this graph. "Va" is the actual wind as taken 
from the soundings and/or streamline analyses. A similar detailed analysis 
of speed differences was made using 210 soundings from areas where 


IGIE exceeded 10 knots. Results are shown in figure 9. 


Wak 











Maps of X for 1000 mb (lower) and 300 mb (upper) are shown in figure 
10. Dashed lines are 5 meter contours, solid lines are 2-kt intervals 
of wind speed, Ix. Figure 11 shows xand Ix for 500 mb and 850 mb. 
These maps were plotted using an FNWF program designed for another 
purpose so that the contour labe!s have no intrinsic meaning. Contour 
interval is 10 meters. Note that the pre-trough region in low levels 
is associated with a x maximum (point of Vx inflow) at upper levels and 
a X minimum (point of Vx outflow) as would be reflected in divergence 
and convergence patterns associated with such height patterns. 

Ix can be seen to be maximum at the 300-mb and 1000-mb levels. 
Note that these maximums are of the same magnitude, 8 to 10 knots, 
but that directions are opposite. As might be expected, speeds are min- 
imum at 500 mb as a reflection of the minimum values of divergence 
known to exist there. 

At 300 mb even maximum values represent a very small percent of 
the wind speed, seldom over 19 percent and generally less than three 
percent of Ve. With the lower wind speeds found at lower levels these 
same magnitudes of Ix are found in some instances to equal 45 percent 
of Ve. Sucha case is depicted along the Gulf Coast in figure 8 (lower). 
It is of further interest that a deepening, fast-moving cold front was 
located in the same area almost directly under the 1000-mb Ix maximum 
center. This front moved rapidly up the East Coast and into the North 


Atlantic. A center of maximum Ix was associated with this front through- 


out the traversal of the United States but speed of movement and the 
short period of accurate data prevented detailed study of this associa=- 
tion. 

Combining Ve and Vx, Vt was formed. Figure 12 (upper) shows the 
effects of this combination at 850 mb; 12 (lower) depicts 1000 mb. 
Solid lines are isotachs of Vt, dashed of Ve. The effect of this 
combination at middle and upper levels is so small as to be undetectable 
in this type analysis and is not shown. The effects at 1000 mb are quite 
evident. Detailed analysis of these proved that here there is a change 
in location of relative differences but no change in magnitude of the 
difference. Terrain and friction effects preclude any meaningful com=- 
parison of the small deviations from Ve to actual wind data. Further 
attempts to contrast Vt with Vg only further emphasizes that there is no 
Significant difference between Vt and Ve. 

Figures 12 and 13 depict the results of using Vg and Ve for advection 
of absolute vorticity. Solid lines are lines of zero advection using Vg, 
DADV; dashed lines are for Ve, EADV. The 300-mb level is omitted but 
showed no marked difference between the two advecting winds as it too 
showed practically congruent patterns. The numerals in the ratios shown 
are located at points of local minimum and maximum value, the upper num- 
eral is for DADV. Units are 10°’ sec ay These values of change of 
vorticity are typical of the period investigated and showed DADV to be 


larger at both extremes by a factor of 5 to 10 units. Computing the pillow 


rs 





and RMSD for all the standard levels gave a pillow O in every case and 
RMSD of 2 to 7 units with largest being at 300 mb. Use of Vt made 

no detectable difference in overall patterns. Some minor increases 

in value in the vicinity of sharp curvature at low levels was noted but 
were of the order of only 10 to 20 percent. Advection of vorticity with 
JX was not computed as its significance has been seen to be small by 
its small effect when added to Ve. However, advection of temperature 


using Vx was computed uSing a program which evaluates the expressions 


J (E,T) (21) 


and 


VeVT . (22) 


One hour change of temperature was computed. Figure 14 shows 1000= 
mb change of temperature with Ve (solid lines), Units are degrees 
celsius lor . Detailed analysis of the grid map shows that advection 
with Vx over the entire map was less than = 0.01 degrees over the United 
States except in the area along the Gulf Coast where values are depict- 
ed with a dashed contour enclosing an area of = 0.2 degrees. Here the 
advection with Vx actually exceeded EADV. This is the area noted earl- 
ier of the deepening low pressure system and front. Relatively large 
values were recorded at the 850-mb and 700-mb levels in this same area. 
Again the rapidity of movement of the front and associated values of ad=- 
vection and short length of the data period prevented detailed investiga- 
tion but suggests the possibility of a meteorologically significant rela- 


tionship. 
14 





9. Conclusions 

1. Where there is curvature in the pressure patterns involved, the non- 
divergent velocity is a better approximation to the wind than the geo- 
strophic value. 

2. No significant change is made by combining the non-divergent and 
divergent winds for the uses investigated. 

3. Arelationship possibly exists between the divergent wind and front- 
al areas that could be utilized in locating such phenomena. This con- | 
clusion is based on very limited data but offers an interesting avenue 


for further investigation. 
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APPENDIX I 


Pillow is defined as 
¥=N=3968 


P= Z_(a-B)n | (I-1) 


n= 





The root-mean-square differnece (RMSD) is determined as follows: 


MSD : > [(a-8)- PIn/x (1-2) 


In equations (I-1) and (I--2), A, and B, represent the point values of 
the fields being compared. In this example the A value is a geostrophic 
value of speed or direction and B is the corresponding value of either non- 


divergent or total wind. 
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